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An exact calculation of the transverse susceptibility for an 
antiferromagnetic Ising A. chain 

Nobutaka Kunisada and Yoshiyuki Fukumoto*) 

Department of Physics, Faculty of Science and Technology, 
Tokyo University of Science, Noda, Chiba 278-8510 

We study the transverse susceptibility of the fully frustrated antiferromagnetic Ising 
Z\-chain, extending Minami's transfer-matrix method for the transverse susceptibility of 
general-type Ising linear-chains [J. Phys. Soc. Jpn. 67 (1998), 2255]. For transverse fields A 
on tip spin sites and i~a on bottom spin sites, we calculate zero-field transverse-susceptibilities 
Xtip = limri.rwo A*S P /A and Xbottom = limri,r 2 -»o Mg ottom /r 2 , where M t ^ p(bottom) denotes 
the magnetization for tip (bottom) spin sites. Both the transverse susceptibilities follow 
Curie's law at low temperatures. We also calculate Xbottom C^i > n )i transverse susceptibility 
of the bottom spin chain under finite tip-spin transverse-fields, to understand the Curie type 
behavior in the zero-field susceptibility. Using the second-order perturbation theory, we 
discuss the A dependence of Xbottom (A) at zero temperature. 

§1. Introduction 

The tranverse-field Ising model has been often used as an effective model to study 
order-disorder transitions in some physical systems. ^ There are many theoretical 
studies for this model. In early time, de Gennes employed this model to understand 
ferro-electric ordering in Pottasium Dihydrogen Phosphate type systems. As for 
the transverse-field Ising chain with only nearest neighbor interactions, magnetic 
properties such as the free energy, magnetization, transverse susceptibility, etc, were 
obtained by Katsura, Pfeuty, and Fisher. 2 )~ 5 ) In those investigations, a mapping to 
a free fermion system via the Jordan- Wigner transformation was done, and thus it 
is difficult to study more general models with S > 1, second neighbor couplings, and 
so on. 

In 1990's, the exact transverse susceptibility for several type Ising chains was 
examined theoretically. Idogaki et al. developed the differential operator method 
to obtain exact transverse susceptibility for a random-bond Ising chain. 6 ) Using 
this differential operator method, Kaneyoshi for the first time obtained the exact 
transverse susceptibility for mixed spin-1/2 and spin-S" (S > 1) Ising chain. 7 ) 

On the other hand, Minami used a transfer-matrix method to calculate the exact 
transverse susceptibility for the Ising linear-chain with arbitrary spin. 8 ) And later, 
he introduced more general Ising- type models with mixed spins and random bonds, 
and calculated the transverse susceptibility for these models by the transfer-matrix 
method. 9 ) Minami's results agree with the previous results obtained by Idogaki et 
al. and Kaneyoshi, and he was also succeeded in getting some original and general 
results. According to Minami's procedure, the calculation of the exact transverse 
susceptibility is reduced to solving the eigenvalue problem of a finite order transfer 
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Fig. 1. Geometry of the antiferromagnetic Ising Z\-chain. The black circles denote tip spins and the 
white circles represent bottom spins. Note that couplings between bottom spins give fourth order 
terms of Fermi operators in the Hamiltonian when we use the Jordan- Wigner transformation. 



matrix. 

In this paper, we consider an exact calculation of the transverse susceptibility 
for an antiferromagnetic Ising Z\-chain under transverse fields i~i on tip spins and 
i~2 on the bottom spin chain (see Fig. [1]). The Hamiltonian is written as 

N-l 

h = jY1 l s 2i4 i+ 2 + 4i+i(4i + 4i+2)\ - AQi - r 2 Q 2 , (i-i) 

i=0 

where are spin-1/2 operators, N denotes the total number of triangles, and the 
periodic boundary condition, si + n = sj, is imposed. The operators Q\ and Q 2 , 
respectively, represent the magnetization operators for the tip and bottom spins 
defined by 

N-l N-l 

Qi = Y.4* + i, Q2=Y.4i- (i-2) 

i=0 i=0 

We now define zero-field susceptibilities as follows: 

Xtip = _ lim MS P /A, bottom = r lim M^ ottom /r 2 , (1-3) 

where M t ^ p = (Qi) and -^bottom = (Qs) are thermal averages of the magnetization 
operators under the Hamiltonian H. 

The present Ising zA-chain, without the transverse fields, is one of the sim- 
plest fully frustrated Ising models with a macroscopic ground state degeneracy, 
and infinitesimal transverse fields determine the ground state uniquely. 10 ^' n ) The 
transverse-field Ising Z\-chain has been known to have a ground-state quantum phase 
transition as a function of the ratio of Zj and rb- 11 ^ We set i~i = r and F 2 = \r. 
The staggered state defined in Fig. [2] (a), which involves spontaneous symmetry 
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Fig. 2. Schematic representations of (a) staggered state and (b) transverse-field paramagnetic state. 
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breaking in the bottom spin chain, is stabilized for A < X c (r), and the transverse- 
field paramagnetic phase defined in Fig. [2] (b) is stabilized for for A > A c (r) . As r 
is decresed, A c (-T) is increased and A c (+0) = 0.984. 

When an infinitesimal transverse field stabilizes the transverse-field paramag- 
netic phase, then M£ ottom > 0. In this case, we expect Xbottom follows the Curie 
law at low temperatures. On the other hand, Mg ottom = in the staggered state, 
and then Xbottom takes a finite value at low temperatures. The main interest in the 
zero-field transverse-susceptibility in the present fully frustrated model is to study 
which behavior is realized in the temperature dependence of Xbottom- As f° r X?ip' 
the low-temperature Curie law is expected because of My > 0. 

In our calculation of the zero-field susceptibilities Xtip an d Xbottom' we use ^ ne 
notation Zj = r and I2 = A-T and write our Hamiltonian as 

H = H -rQ {Q = Qi + AQ 2 ), (1-4) 

where 

N-l 

h = jY1 \4A l+ 2 + 4i+Mi + 4i+2)\ (i-5) 

i=0 

is the unperturbed Hamiltonian. Because lim/^o^ = Hq has a diagonal form, 
the calculation of transverse susceptibilities can be made by a direct extension of 
Minami's method. 9 ) In §2 we show both Xti p an d Xbottom diverge for T — > 0. 

The divergence of Xtip is reasonable as mentioned previously, but the divergence 
of Xbottom I s n °t trivial. In order to get further insight into the low temperature 
behavior of Xbottom' we k ee P A finite and study the susceptibility Xbottom (A)- ^ n 
this case, we write our Hamiltonian as 

H = H - r 2 Q 2 , (1-6) 

where 

N-l 

n = J i s 2i4i+2 + 4 i+ i(4i + 4i+2)} - AQi- (i-7) 

i=0 

It should be noted that TCq has off-diagonal terms, so we have to start with the 
diagonalization of 71q. This procedure is explained in §3. As a result, we find that 
the above-mentioned divergence is suppressed by i~i > 0. 

Needles to say, exact solutions of one-dimensional quantum spin systems have 
provided testing grounds for the development of numerical calculation methods. 12 )' 13 ) 
A typical example is the quantum XY model with nearest neighbor inter actions. 2 ) 
In most of such systems, however, the structure of energy spectrum is so simple that 
a numerical calculation method reproduces the exact results rather easily. Therefore, 
it is desired to find exact solutions of frustrated systems with complex energy spectra. 
Our exact transverse susceptibility of the fully frustrated Ising Z\-chain may be useful 
from this point of view. 
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The present paper is organized as follows. In §2, we calculate Xti P an d Xbottom 
by using Minami's transfer matrix method. In §3, the calculation of Xbottom (Ti ) * s 
presented, where we use an unitary transformation introduced by one of the present 
authors. 11 ^ In §4, we discuss the behavior of the zero-temperature susceptibility as 
a function of r% on the basis of the second-order perturbation theory. The results 
obtained in this paper are summarized in §5. 

§2. Exact calculation of zero-field susceptibilities Xtip an d Xbottom 

2.1. Formulation of the transverse susceptibility 

The transverse susceptibility \ x f° r the system described by the Hamiltonian 
(jl-4p is defined as follows: 

X 



dr Tre-P H 



d_ 

dr 



r=o 



(Q)_ + rf P dt(e tH °Qe- tH °Q) 



p +°( r 



i + rf p dt(Q) 



r=o 



o 



dt(Q(t)Q) Q , (2-1) 



where f3 = 1/T is the inverse of temperature, (• • -) = Tr • • • e P Ho /Tre P H ° : and 
Q(t) = e tH °Qe~ tH ° . In Eq. (|2-ip . we have used the well known formula, 



e -0H ^ e -/3H 



jl + T J dtQ(t) + ---\, (2-2) 



and the fact of (Q) = 0. Then the transverse susceptibility for the antiferromagnetic 
Ising Z\-chain is given by 

X X = Xtip + ^Xbottom ( 2 ' 3 ) 

with 

Xfip= f dt(Q 1 (t)Q 1 ) = N [ dt(sf(t)sf) o , (2-4) 
Jo Jo 

Xbottom = [ dt (Q 2 (t)Q 2 ) = N [ dt (s x Q (t)s%) , (2-5) 



where we have used the translational invariance. 
2.2. Transfer matrix 

Now we introduce the transfer matrix for Eq. (j 1 • 5 p . Let \l) be the eigenstate 
of s z with the eigenvalue di (d± = —d 2 = 1/2). Then the transfer matrix V can be 
written as follows: 

2 

V m ,n = £ V Un {yUn = e~^ d ^ +d ^ +d ^) . (2-6) 
1=1 
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In a matrix form, V is given by 

_ /2e-W4 cos h(/3 J/2) 2e^/ 4 \ 

^"V 2e' 3J / 4 2e-/ 3J / 4 cosh(/3J/2) y ' • ^ 

The maximum eig envalue of V is Ai = e" 3/3J / 4 (l + 3c 13 - 1 ), and the minimum one is 
A2 = e _3/3J//4 (l — e@ J ). Eigenvectors of V are obtained as follows: 

Then, we can define an unitary matrix 

U=( U11 Ul2 )=^(] \)=rf, (2-9) 

\U 2 l «22/ v/2\l "V 

which diagonalizes V as A = WVU, where A is a diagonal matrix whose elements 
are (A)^ = Aj<%. 

For an integer L = O(N), we calculate V L as follows: 

V L = U^A L U=( UU U21 )( X J Ul2 )^^-(] ]) (2-10) 

\U21 U 12 J \0 \%)\U21 U22) 2 V 1 V 
in the thermodynamic limit, N — > 00. The partition function is 

Z = Tte- f,Ho = TrV N = TrA N ^ Af. (2-11) 

2.3. Exact expressions for the zero-field susceptibilities 

To obtain x£ p and Xbottom' we nee d to calculate correlation functions, (sq(*) s o)o 
and <sf(*)af) . Defining /»§(») = J[s z 2i s z 2i+2 + s z 2i+1 (s z 2i + s z 2i+2 )}, we rewrite (sf (t)sf} 
as follows: 



^o-^^Ki 1 5* 1 0l 2 W£ -1 . ( 2 - 12 ) 



u p,q=ll,l' = l 

where |{/}) = |Z ,-" >fejv-i), = E? =i " " " £L-i=i> and ^(0 = J ( d P + 

Similarly, we obtain 



(s°(t)$) = ±- V ({I }| e *^(- 1 )+^(o)] s g e -*W(- 1 )+hS(o)] s g e -^o | {Z }) 

= 4- E E et[Ws(/) " Ws(/ ' )] K^ I s * 1 OlXx^" 2 (2-i3) 



p,q,r,s=ll,l'=l 
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with £ pqrs (l) = J(d p + d q + d r + d s )d[. The transverse susceptibilities per site, 
X0,tip = lim A^°c Xtip/ N and A0,bottom = limjv^oo x£ ottom /iV, turn out to be 



Xo,tip 



4A 



-/SJdpdq 



(2-14) 



where A£ pq (l) = £ pq (l) — £ pq (l) with Z = 3 — Z, and 



1 



- 0Spqrs {I) 



Xo,bottom - \2 ztw /ir m 



j3Spqrs {]>) 



-0J(d p d q +d r d s ) 



, (2-15) 



ASpqy* s ^0 



where /\£ p q TS {t) — £ p qrs(J^) £pqrsi] J }' 

We carry out summations in Eqs. (|2-14p and (|2-15p . and then obtain the exact 
transverse susceptibility of tip spins and the bottom spin chain as follows: 



Jxo,ti P - 2Ai 

Jy* -J- 

J AO, bottom ~~ ^2 



e /3J/4 _ e -3/3J/4 + p Je PJ/± 



-f>J/2 _ e -WJ/2 + 2/3J / -/3J/2 + 2e /3J/2\ 



(2-16) 



, (2-17) 



which are plotted in Fig. [3J When temperature decreases, Xotip increases and di- 
verges as expected. The low-temperature asymptotic form of Xo tip ^ s g iVen by 



Xo.tip 



1 t 1 
6T + 6J' 



(2-18) 



The Curie constant is 1/6, which is comparable to the value in the high-temperature 
limit, 1/4. (See the inset in Fig. [3j where Xotip^ ^ s shown.) We also find that 




Fig. 3. The exact transverse susceptibility for the antiferromagnetic Ising Zi-chain. Equations 
(|2-16[) and (|2-17|) are plotted as a function of T. 
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Xo.bottom diverges as T -> 0. Looking into the expression of X§,bottom in Ec l- uM> 
then we notice that it is not affected by the transverse field T\ at all. Therefore, 
the corresponding ground state in the calculation of \q. bottom * s the- transverse-field 
paramagnetic state shown in Fig.[2](b), and thus the low-temperature Curie behavior 
appears in Xq bottom due to a finite bottom-spin magnetization along the x direction 
induced by the infinitely small r 2 at zero temperature. The asymptotic form of 

Xo,bottom * S 

Xo,bottom — ^ + 7~j- (2*19) 

The Curie constant 1/9 is somewhat smaller than that of Xotip - ^ ne t em P era t ure 
dependence of %q. bottom-^ 1S shown in the inset in Fig. [3] 

In the above susceptibility calculations, we have set 7i, r 2 ~~ * at the same 
time. If we keep A finite and calculate Xo bottom (A) 5 then the corresponding ground 
state turns to the staggered state shown in Fig. [2] (a) and different low-temperature 
behavior is expected. In order to study this issue, we carry out an exact calculation 
of Xo. bottom C^i) m ^he following section. 

§3. Exact calculation of transverse susceptibility Xbottom(^i) 

3.1. Unitary transformation ofTCo 

In this sectoin, we consider the Hamiltonian in Eq. (jl-6j) for finite values of 
r%. In order to carry out the exact calculation of the transverse susceptibility, we 
have to diagonalize TCo . The diagonalization is achieved by introducing the following 
operator: 

2 

Ri= J2 Pi{d*,dpy 2e ^ d ^+\ (3-1) 

a,/3=l 

where 

Pi(d a , dp) = \(l + 4d a a| j )(l + M p s z 2l+2 ) (3-2) 

is the projection operator for the two bottom spins on the ith triangle. The angle 6 
is defined by 

e(di,dx) = -9(d 2 ,d 2 ) =tan- 1 ^^, e(d 1 ,d 2 ) = -6(d 2 ,d 1 ) = j, (3-3) 

1 1 4 



where Q = \J J 2 + if*. This angle represents rotation angle of the tip spin depending 
on the bottom spin state. n ) (See Appendix for the derivation of Ri.) Then, the 
unitary transformation R diagonalizing Eq. fjl • Tj) is 



N-l 



R=]jRi. (3-4) 

i=0 

Carrying out the unitary transformation by R, we obtain the transformed Hamilto- 
nian H as follows: 

H = R^HR = H Q — r 2 Q 2 , (3-5) 
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N-l 



N-l 



HO = rfn R = J £ 4i4i+2 + E ( J + S 2i + J-4 t+ 2) 4i+l, 



(3-6) 



i=0 
N-l 



i=0 



Q 2 = rfQ 2 R = £ [J 0S |. + hslAi-Ai+l + I+4i4i+l - I-4i4i-l] . (3-7) 



i=0 



where the coefficients are defined by 



T -O+T T - J T - 2J T - ° ±Fl 



(3-8) 



3.2. Transfer matrix 

Let us introduce the transfer matrix V for T^o- As the same way in §2, it is given 

by 

2 

V — V^V' fv' — -/3[Jdmd n +(J + d m +J-d n )di]\ (1.Q\ 

Vm,n — y m,l,n \ y m,l,n — e J 3 1° y 7 



or, in a matrix form, 



V 



2e" /3J / 4 cosh 2e^ J / 4 cosh /?A 

2e /3J/4 cosh p Fl 2e - /3J / 4 cosh y ' 



(3-10) 



The maximam and minimum eigenvalues of V are, respectively, Ai = 2e~^/ 4 cosh 
2e^ J / 4 cosh j3F\ and A 2 = 2e~ /3J ^ 4: cosh j3 17 — 2e^ //4 cosh /3-Ti- The corresponding 
eigenvectors are those in Eq. (|2-8p . The partition function is 



Z n = Tr e 



(3-11) 



3.3. Exact expression for Xbottom(-^i) 

The transverse susceptibility for H in Eq. (|3-5p can be written as 

XbottanCA) = ^ dt(Q 2 (t)Q 2 ) Q , 

where 



(3-12) 



Q 2 (t)Q 2 ) o = v 



4 <*?(*K>o + A 2 <[s§«](^§«} 
+4 <M(W> + /2 - 



(3-13) 



In the above equations, we have redefined as (• • -) = Tr ■ ■ ■ e ^ w oy"pr e Q 2 (t) 
e mo Q 2 e' m ° and so on. In the thermodynamic limit N — > 00, we obtain 



xS,bottom(A) = hm XbottamCAj/JV 

iv^oo 



lim 



• (3-14) 
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It is convenient to define h^i) = Js^s^,^ + (J+$2i + ^- s 2i+2) s 2i+i- Then, the 
correlation functions are written as follows: 

{s*(tX) = i J2 d 1 >l e*^(- 1 )+^(0)] s g e -*[^(-i)+^(o)] s g e -^o | { j }) 

= j- E E et[4,rsy ° )_49rs( " o)] Ifa I I Zo>l 2 v;^x, r Xp~ 2 ( 3 - 15 ) 

P,q,r,sl 0l l' Q 

with £ pqrs (lo) = J(d p + d s )d io + J+d r d k) + J~d q di , 

(*)*S«> 

= i e *[4«(fe.»-i,»i)-4i(»i.»'-i,ii)] |(/ | s g | /^)| 2 

p,q lo,h,l-i 
V V V 

x K'-i I * y -i I l '-i)\ 2 \(h I s{ i zoiX^x^Xr 2 (3-16) 

with Spq(l , 1-iJi) = J(d p + d q )di + J + (di Q d h + + J^(d h d p + d^di^), 

= _L £ ({ z }| e t[^(-i)+^(o)] a v a « e -t[&S(-i)+hS(o)] a v a v e -/9Ho |{ z }) 



iV-2 



I' V 

(3-17) 

with £+ qr (lo, h) = J(d p + d r )di + J+di d h + J-(d h d p + dj d g ), and 

= J_ V ({/ }| e 4^(-i)+^(o)]^^ ie -4^o(-i)+^o(o)]^^ ie -^o | { ^ }) 
F" E ^ et[ ^ r(, ° ,Ll) ~^ r(,&,, - l)] |<^0 | si | ^}| 2 |<Z_! I si, I Z^ x > | 2 

1 0>'-1 

xV; Ll ,i V^ ?jr V^- 2 (3-18) 



^ _ . __ 

p,q,rl ,l_ l 



with £ pqr (lo, l-i) = J(d p + d r )di + J + {dpdi_ 1 +di d q ) + J-d^di^. 

We calculate summations in Eqs. (j3-15[) - f)3-18n . and substitute those results into 
Eq. (|3-14p . Carrying out the integration over t, we obtain our final expression: 

Xo,bottom(A) = IqCxx + lfC X yy + I+Cy y + I_C yy . (3T9) 
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Here, C xx is given by 

r P 

C.r 



lim / dt(s x (t)s x ) 

N— >oo ' 







4A? 



-j3{J+d v dq+J-d r d a ) 



1 p,q,r,s E pqrs (l) £pgrs(2) 



2A< 



e- pil (fj+n + fj-a) + e-^Vj-A + /j+a 
+2e -/3(J-«+2r 1 )/2 /fl cosh + 2e -/3(J + r l)/2 /a cogh 



(3-20) 



where A' x = [1 + e"^ 1 + 2e" /3 ( J+ri )/ 2 cosh/?J?] 2 and = (1 - e~ Px )/x. Similarly, 



C X yy denotes 



C xyy = lim f dt{[sts y _ lS lm4s y -A), 



1 



E E 



/3^pg(io,i-l,^l) 



64 ^l lof^^Pli^J-lih) ~ Epq(lo,l-l,h) 
_ C X x 

~ 16"' 

where we have used the notation Fj = 3 — Zj. Finally, C~ denote 



(3-21) 



N— >oo 
1 



C±= hm / d^4iK*K4i> 



16A 
1 

8A[ 



EE 



P[£pqr(loJ) + J±dqd r 



+ 



/3 



p,q,r l ,l 
■± 
P<1 



A£p qr (lo,l) 32A 2 



A£pq r ^0 



EE-" 

p,g,r lo,l 

Af^ — 



+2e -/ 3 (j-i2 + 2r 1± r l) /2 ; + 2(3e -p(j+2n)/2 cogh gC£±A) 



, (3-22) 



where A£^ qr (lo, I) = £^ qr (lo, I) — £^ r (Zo, 0- If we se * A — > in our final expression, 
we obtain 



lim J X X 

.bottom (A) 

1 1— >u 



5 _ 4e" /3J - e" 2/3J + f3J{2 + e" 
2(3 + e-' 3 - 7 ) 2 



(3-23) 



which is the same as Xq bottom obtained in the previous section. 

In Fig. H] (a), we show the temperature dependence of X§ bottom (-^i ) f° r = 
0,0.2,1,2,10 and oo in the unit of J = 1. It is found that %g bottom (A > 0) has 
a plateau at low temperatures instead of the Curie law of Xq bottom (A = The 
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appearance of such plateau structure is naturally expected if the staggered ground 
state is realized. The most transparent situation in this context is the large A limit. 
Setting A — » oo in our expression of Xo bottom C^i)> we obtain 



lim Jxg 



0, bottom 



/3J 



cosh^(/3J/4) 



+ tanh(/3J/4) 



(3-24) 



which is nothing but the transverse-field susceptibility of the antiferromagnetic Ising 
chain with nearest neighbor interactions. In the limit A — ► oo, the tip spins polarize 
to the x direction completely and the remaining degrees of freedom are just bot- 
tom spins, which follow the nearest-neighbor Ising model. For the nearest-neighbor 
Ising chain, it has been well known that such a plateau structure appears in the 
temperature dependence of the transverse susceptibility. 9 ' 

It is interesting to study the zero-temperature value of x% bottom C^i)- Assuming 
A 7^ 0, we take the limit of T — > in Eq. (|3T9p . Then we get 



J™ Xo,bottom(A, 



1/4 1 1 
8 1 J + A + J + A 



(3-25) 



which is plotted as a function of A in Fig.[J](b). When A — > 0, the zero-temperature 
susceptibility diverges like 1/8-Ti. On the other hand, it approaches to 1/2 J for large 
A- 



(b) 5 : 




Fig. 4. (a) Temperature dependence of Xo,bottom(A) for several values of A, and (b) zero- 
temperature value of Xo bottom (A) as a function of 7\, where we set J = 1. 



§4. Zero-temperature transverse susceptibility: perturbation study 

In the previous section we have obtained the exact expression of the transverse 
susceptibility Xq bottom (A) an d found the low-temperature plateau structure for A > 
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0. In order to get deeper understanding to this behavior, we here use the second-order 
perturbation theory to calculate the zero-temperature transverse susceptibility. 

The macroscopic degeneracy of the Ising Z\-chain is removed by and we 
obtain the unique ground state </>g) of as shown in Fig. [5l which is, of course, 
equivalent to the staggered state. The ground-state energy is 



N 



(J + J+- J- 



(4-1) 



The perturbation operator Q 2 in Eq. (|3-7p is applied to this ground state, then we 
obtain four types of intermediate states, \ipi) (i = 1 ~ 4), with the energy of e,;, 
which are defined in Fig. The energy denominators, AEi = e% — eg, are 



AEi = AE 2 = J + 



~2 



AEi = J 



J_ 

T 



AEa = J + 



(4-2) 



2 2 ' 2 ' 2 

Thus, ground-state energy Eg(r 2 ) perturbed by the transverse field r 2 is given by 



Eg(r 2 



J 2 



+ 



8(J + A) 16(J-fl + A 
The zero-temperature susceptibility is obtained by 

d 2 



+ 



(rt-A 



16(J + ^ + A) 



.(4-3) 



Xo , bottom (A 



T=0 



or 2 



Eg 
N 



(4-4) 

which leads to Eq. (|3 25|) . Let us consider the limit i~i — > 0, where Xo bottom (A) — 
I/8/1. Noting Z\i?3 — 0(/i), we find that the divergent behavior for Fi stems 
from the intermediate state (^3), which is a member of the degenerate ground-state 
manifold of the Ising Z\-chain with no transverse fields. (Note that the number of 
antiferromagnetic bonds in 1 9^3) is the same as </?g)-) In other words, introduction of 
finite tip-spin field i~i opens an energy gap between \tpg) and \(fs) and thus prevents 
the divergence of Xobottom(A) at low temperatures T <~ r%. 



\tpg) , eg (ground state) 
2i+\i 
J+f* J. 



\<Pl) - 



kg). 



\f2) ■ 



% s 2i+l S 2i-l |Vg; . f 2 




AA AA 



M = 4 s Ii+ikg)> £3 



A\ AA 



Fig. 5. Spin configurations of the ground state \<pg) with energy eg and intermediate states \<pi) with 
the energy (j = 1 ~ 4) in the second-order perturbation theory. The strength of the exchanges 
is J + > J > J_. The black circles represent spins flipped by a perturbation operator. The black 
and gray bonds denote antiferromagnetic (AF) and ferromagnetic (F) bonds, respectively. 



An exact calculation of the transverse susceptibility 



13 



§5. Summary 



In this paper, we have studied the transverse susceptibility of the antiferro- 
magnetic Ising Z\-chain. The spin-1/2 antiferromagnetic Ising Z\-chain has next- 
nearest-neighbor interactions, so it can not be mapped to a solvable model by the 
Jordan- Wigner transformation. We have employed Minami's method to carry out 
an exact calculation of the transverse susceptibility of the Ising Z\-chain. 

First, we have considered the system on which infinitesimal fields i~i on tip spins 
and i~2 on bottom spins are applied. Using the transfer-matrix method, we have 
obtained the transverse susceptibilities, Xotip f° r A an d Xo bottom f° r A- We have 
found that both xg tip and Xo. bottom shows Curie type behavior at low temperatures. 

Next, we have extended our calculation to the the antiferromagnetic Ising A- 
chain under a finite field F\ on tip spins and an infinitesimal field i~2 on the bottom 
spin chain. We have diagonalized the Hamiltonian by an unitary operator R. In- 
troducing the transfer matrix for this transformed Hamiltonian Hq, we have carried 
out the exact calculation of the transverse susceptibility for i~2 which is denoted by 
Xo bottom (A)- We have found that finite T\ changes the low-temperature behavior 
and Xo. bottom (A) shows the low-temperature plateau. The role of i~i has been also 
discussed by using the second-order perturbation theory, by which the spin configu- 
rations giving rise to the low-temperature Curie law have been identified. 
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Appendix 

In this Appendix, we describe how to introduce the operator Ri in Eq. (|3T 
Let us rewrite the Hamiltonian Eq. fjl • Tj) as follows: 

N-l N-l 



Ho = JY, S 2i4i+2 + S 2i+ 2 ), (A-l) 



i=0 i=0 

where 



hi{S2i, s 2i+2) — J s 2i+l\ s 2i + s 2i+2> ^l s 2i+l 

_ip(4 + 4 +2 ) -A 
" 2 V -A + 4 i+ 2, 



(A-2) 



In the second line of Eq. (|A-2p . we have used the matrix representations of and 

S 2i+1- 

When s z 2i = s|j +2 = 1/2, the eigenvalue equation for the matrix in Eq. (IA-2P is 
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Eigenvalues are calculated as 



e = ±\^P + r^±^, (A-4) 



and eigenvectors are 



a ^ 1 f r t \ ( cose ^ fore = ± |_ (A . 5) 



v yjr( + j) 2 vf(^ - J )J v^ sin 



From Eq. (|A-5p . we can obtain an unitary matrix Ri(d\,di) which diagonalizes 

t-> / r i \ ( cos 6 sin#\ „ v . _ ^f) q y ,. 

B 4 {dx,d 1 )= [ sin0 cos J = cos + *2 S | +1 sin 9 = e 2l6s ^ . (A-6) 

In the same way, we can also obtain the unitary matrix i?i(d2,^2) for sfi = $21+2 = 
— 1/2. The eigenvalue equation 

leads to the following eigenvalues and eigenvectors: 

a = —sin 9, b = cos 9 for e = — , (A-8a) 

a = cos#, b = sva.9 for e = — — . (A-8b) 
Thus, the unitary matrix ^,(^2,^2) diagonalizing ^(^2,(^2) is 

/ , ,\ /cos 6* — sin#\ — 2?(9s y /a ^\ 

=U* cos^J =e 28+1 - ^ 

For the case of sfj = — sfj , 2 , the eigenvalue equation 



gives 



-6=i= fore = | (A-lla) 



a = b = for e = - — . (A-llb) 
v 2 2 



Then Ri(d\,d2) and Ri(d2,d\) are, respectively, given by 



1 / 1 1 



i^di.da) = ^ 1 ) = e ^ 2i+1 > ( A42 ) 
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R i (d 2 ,d 1 ) = ^(] i) = e-^w. (A-13) 

From Eqs. (|A-6p . (|A-9|) . (|A-12p and (jA43j), we can define Ri using the projection 
operator (|3-2p as follows: 

2 2 

Ri = ^ y £ J Pi{d a ,dp)R i {d ol ,dp), (A-14) 

a=l f3=l 

which is Eq. (13-11) . 
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